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Abstract. Originating in different depths of the very extended atmospheres of AGB stars, various molecular spectral lines 
observable in the near-infrared show diverse behaviours and can be used to probe atmospheric dynamics throughout the outer 
layers of these pulsating red giants. In Nowotny et al. <20051 Paper I) time series of synthetic high-resolution spectra were 
presented, computed from a dynamic model atmosphere for a typical carbon-rich Mira. In this work, line profile shapes, their 
variations during the lightcycle and radial velocities derived from wavelength shifts are analyzed and compared with results 
from observed FTS spectra of the C-rich Mira S Cep and other Miras. It is found that the global velocity structure of the model 
is in qualitative agreement with observations. Radial velocities of molecular lines sampling different layers behave comparably, 
although some differences are apparant concerning absolute values. A correction factor of p » 1.36 between measured RVs and 
actual gas velocities is derived for CO Av=3 lines. It is shown that dynamic model atmospheres are capable of reproducing 
Mira spectra without introducing an additional "static layer” proposed by several authors. 
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1. Introduction 

High-resolution (IR) spectroscopy has been successfully used 
to study the kinematics within all regions of the extended at¬ 
mospheres of AGB stars, as summarised in detail in Nowotny 
et al. J2005I in the following designated Paper I). The outer 
layers of these red giants are heavily influenced by pulsation 
of the stellar interior and by the development of a dust-driven 
stellar wind. From observations we know that different molec¬ 
ular lines sample layers of various depths and therefore show 
typical behaviours depending on the movement of gas at the 
corresponding optical depth. 

Based on a dynamic atmospheric model, synthetic spectra 
were calculated by applying a radiative transfer that includes 
the influence of velocities. These were presented in Paper I, 
demonstrating that the models can qualitatively reproduce dy¬ 
namics throughout the different regions of the atmosphere. 
Here, the results of our modelling (line profile shapes and their 
evolution as well as measured radial velocities) will be quan¬ 
titatively compared to observations, especially to the ones of 
S Cep. 


Send offprint requests to: W. Nowotny, 
nowotny @ as tro. uni vie. ac. at 


2. Observational results of S Cep 

S Cep is the only carbon-rich Mira for which extensive time 
series IR spectroscopy has been obtained so far. As this pa¬ 
per deals with C-rich model atmospheres, S Cep was therefore 
chosen as the reference object to define the model parameters 
for our calculations. Observational results were published by 
Hinkle & Barnbaum (I1996I in the following HB96). S Cep is 
the brightest C-rich Mira in the northern sky and is believed 
not to be significantly obscured by dust (HB96). 

HB96 obtained high-resolution spectra in the range of 1.6- 
2.5;um by using the FTS spectrograph at the 4m Mayall tele¬ 
scope at Kitt Peak. The spectral resolution was set to R«i70000 
and the S/N ratios are of the order of 50. Details can be found in 
their Table 1. In addition to the published K-band radial veloc¬ 
ities (RV), the spectra were kindly provided to us by K. Hinkle 
for the present comparison with our models. All spectra were 
corrected for earth motion and the stellar velocity' was sub¬ 
tracted. 

Absorption lines of the CN molecule have been used in the 
past to study velocities in C star atmospheres (see references in 
Sect. 2.3 of Paper I) as they are prominent in spectra of carbon 
stars. As the crowding of CN lines is very high, HB96 used 
cross-correlation of a region of 50cm“’ width in the K-band 

' The center of mass radial velocity (CMRV) was derived from pure 
rotational molecular (CO, SiO, CS, HCN) emission lines in spectra 
observed with radio telescopes, see TableU] 
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to measure RVs, instead of line shifts of single lines. Also sam¬ 
pling deep photospheric layers, RV-curves from CN lines in the 
near IR (Fig. 1 of HB96) show a similar behaviour to the ones 
for CO Av=3 (Paper I). This was also found by Alvarez et al. 
( 12001> from cross-correlation of visual spectra of S Cep (their 
Fig. 3). 

Figure[2shows observed RVs taken from the literature (de¬ 
tailed explanations in HB96) as a function of visual phase (py. 
All measurements from different light cycles are combined and 
then repeated for better illustration. While CN Av=-2 lines of 
the red system show the typical discontinuous RV-curve, CO 
Av=2 low-excitation lines vary less and not regularly. The KI 
line, observable in the visual at a;0.77yum, probes outflowing 
regions further out in the atmosphere (as CO Av=l lines) and 
shows RVs comparable to radio measurements of the circum- 
stellar expansion (TableQ. 

It is remarkable that the RV-curve of CN lines is asym¬ 
metric relative to the CMRV and extends to more negative ve¬ 
locities, while all observed CO Av=3 RV-curves (e.g. Fig.l of 
Lebzelter & Hinkle IT002t appear shifted to more positive val¬ 
ues. 


Table 1. Observed and derived properties of the carbon-rich 
Mira S Cep, taken from the literature: 1 - from P-L-relation 
of Groenewegen & Whitelock (I1996II . 2 - and mboi from 
Bergeat et al. (1200 H -i- Hipparcos distance, 3 - mboi by interpo¬ 
lation/integration of VJHKLM and IRAS magnitudes (as de¬ 
scribed in Kerschbaum & Hron[T996j -i- Hipparcos distance, 4 
- GCVS catalogue, Kholopov et al. lll988t . 5 - Neugebauer & 
Leighton lll969t . 6 - Whitelock et al. ( I2000t . 7 - Hipparcos cat¬ 
alogue, 8 - van Belle et al. d 19971 . 9 - Eglitis & Eglite d 19971 . 
10 - Olofsson et al. ( fT^ . 11 - Bujarrabal et al. (dHUl, 12 - 
Loup et al. dl9931 . 13 - Olofsson et al. dl9931 . 14 - Neri et al. 
lfT9^ . 15 - Gautschy-Loidl et al. (|2004j. 
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Fig. 1. Observed radial velocities derived from different lines 
in high-resolution spectra of the carbon-rich Mira S Cep, com¬ 
piled from HB96, Barnbaum d 19921 and Olofsson et al. dl9931 . 
Eor S Cep. v —v/j^//oceHf"("26.65, Vheiiocent—^LSR 11.65 [km s ] 


3. Stellar parameters: real stars vs. models 

Tabled lists stellar parameters of S Cep derived from obser¬ 
vations, which can be compared with the ones of the model 
atmosphere used in Table 1 of Paper I. Relating these sets of 
parameters is not straightforward. In this respect, the dynamic 
atmospheric model used here has to be regarded as having jus¬ 
tifiable and physically reasonable properties rather than being a 
dedicated ”ht” to the C-rich Mira S Cep. Gautschy-Loidl et al. 
( EUn?b were able to reproduce low-resolution spectra (ISO) 
rather well on the basis of this model and in this paper we fur¬ 
ther analyse the predictions of this model. 

3.1. Period, luminosity, amplitude, mass 

Only very few parameters are clearly derivable from observa¬ 
tions and can directly be adopted for the modelling, like the 
pulsational period P (490“^ was used for the modelling). The 
luminosity L can be constrained to some extent (Hipparcos dis¬ 
tances). In this respect, the chosen mean luminosity of 10"^ L© 
for the model is plausible (although probably at the high- 
luminosity end of the range). The amplitude Amboi of 0.86™ 
represents a realistic value comparable to observations of Mira 
variables (e.g. Whitelock et al. Eunnt . Only rough estimates 
can be given for stellar masses of galactic held stars such as 
SCep. The chosen AGB mass of M»=lMo is reasonable, al¬ 
though some statistical arguments (Barnbaum et al. ll991t point 
towards main sequence masses in the range of 2.5 ^Mq for 
certain C stars, including S Cep. In Sect. l6.1l we also discuss a 
model with 1.5 M©. 
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3.2. Effective temperature, surface gravity 

Relating the classical parameters Teff and log g of observations 
to our models is hampered by two fundamental difficulties. 

The first one is inherent to the dynamic models themselves. 
As described in detail in Hofner et al. (I2003> the model calcula¬ 
tion starts with a hydrostatic initial model. Being very similar to 
other model atmospheres for cool giants (as e.g. the ones com¬ 
puted with the MARCS code), they are characterised by the 
following fundamental parameters: luminosity L,, mass M,, ef¬ 
fective temperature T, and the elemental abundances (the C/O 
ratio is especially important, all other abundances are assumed 
to be solar). These can be directly linked to the correspond¬ 
ing standard parameters of stellar atmospheres (Teff, L, log g, 
chemical composition). Then stellar pulsation comes into op¬ 
eration, simulated by a variable inner boundary (piston), and 
dynamical effects (extension, shock fronts, dust-driven winds) 
occur. The atmospheric structure is thus fundamentally differ¬ 
ent to the initial model at any point in time (compare Fig. 9 in 
Paper I or Figs. 3 and 4 in Hofner et al. 120031 . Not only is the 
spectral appearance thereby changed ( Gautschv-Loidl I200 11 . 
but also the structures can no longer be characterised by the 
parameters T, and R* (or log g). In a more general context 
Baschek et al. flWfl assert that the classical term of ’’effec¬ 
tive temperatures” is highly questionable for Mira stars. This 
topic was also discussed by Bessell et al. ifT^ . 

The second problem is related to deriving values for Teff 
and log g from observations of pulsating giants. The parameters 
available in the literature rely on interpreting observational data 
with the help of models (fitting synthetic spectra to observed 
ones). This process always implies uncertainties. This is es¬ 
pecially so, as all studies in the past utilised hydrostatic model 
atmospheres. These may resemble AGB stars with only modest 
variability and mass loss quite well (e.g. Jprgensen et al. l2000t . 
but obviously are not an adequate approach for evolved red gi¬ 
ants with more pronounced pulsations and strong outflows (as 
outlined e.g. in Hron et al. l2()()2l or Gautschy-Loidl et al. l2()()4t . 

Thus, in the end a direct comparison of Teff and log g, de¬ 
rived ’’inappropriately” from observations, with initial model 
parameters exhibiting a very indirect link to the actual dynamic 
model structure is somewhat misleading. 

In the study presented here and in Paper I we use a dynam¬ 
ical model atmosphere possessing a T* of 2600 K, which ap¬ 
pears somewhat higher than the values for Tefr derived from 
observations of S Cep (Table[3- But - this being only a param¬ 
eter of the hydrostatic initial model - it is not a colour tempera¬ 
ture of the dynamic model and therefore not directly compara¬ 
ble. Despite some discrepancies, Gautschy-Loidl et al. (I2004> 
showed that synthetic spectra, based on the same model, resem¬ 
ble observed low-resolution spectra of S Cep in the IR rather 
well. 

3.3. Eiementai abundances, C/O ratio 

For C/O only the second problem discussed in Sect. l3.2l has to 
be considered. Gautschy-Loidl et al. J2004t used the C3 fea¬ 
ture at 5.1 yum as an indicator of C/O. The estimated value of 
C/0=L4 (from the modelling and also from their comparison 


of observed spectra of several stars) appears relatively plausi¬ 
ble. Keeping in mind the limited accuracy (around 0.2-0.3) in 
determining this quantity, the ratios given in Table[2seem to be 
relatively consistent. 

3.4. Consequences 

In spite of the difficulties in relating (Teff, log g) to the funda¬ 
mental parameters of dynamic models, the derivation of these 
parameters is definitely desirable. A major step in this direc¬ 
tion has been made by Gautschy-Loidl et al. J2004t . who at¬ 
tempted to connect dynamical models and observations of real 
AGB stars for different phases by fitting spectra (e.g. ISO) and 
photometric indices. Starting from a grid of models they arrive 
at a ’’best fit model” for a given star by (fine-)tuning the model 
parameters. Although they list the resulting fundamental pa¬ 
rameters [T,, Lt, M,, C/O) of the corresponding hydrostatic 
initial model, Gautschy-Loidl et al. J2004t emphasise that this 
’’fitting” cannot (yet) be used to determine the stellar parame¬ 
ters for a given star. 

To proceed in this matter, further observational constraints 
have to be taken into account, for example detailed line profile 
studies (this paper) and intensity profiles/visibilities (e.g. Hron 
et al. l2003t . Due to the complex structures of Mira atmospheres 
only such a multi-dimensional approach will lead to the desired 
results in the near future. 


4. Comparison of spectra and radial velocities 

We compare both line profiles and radial velocitiy variations 
in observed and synthetic spectra. The plotted observed spec¬ 
tra are normalised with respect to a nearby pseudo-continuum 
(for a discussion of FTS data reduction see Hinkle et al. ll995> . 
Synthetic spectra, on the other hand, are normalised relative 
to a computation with only the continuous opacity taken into 
account. In addition, a wavelength-independent value of opac¬ 
ity is added for the contribution of dust and molecules affect¬ 
ing the spectra pseudo-continuously (cf. Sect. 5.1. in Paper I). 
A comparison of line depths of observations and calculations 
may therefore suffer from some uncertainties. 

Synthetic radial velocity variations were derived from the 
shifts in wavelength of the resulting line profiles. This was done 
by: (i) direct measurements of the deepest point of the line, (ii) 
fitting splines through single profiles and measuring the mini¬ 
mum, (iii) building the first moment of the flux distibution of 
the line or (iv) cross-correlation of a whole spectral region us¬ 
ing the IRAF task fxcor. HB96 used a spectrum of BL Ori 
(C star with only small variations) as a template for the cross¬ 
correlation. Here, a spectrum for maximum light was used, 
computed without taking velocities into account. Synthetic ve¬ 
locities are then compared to observational results. 

For the characterisation of radial velocities we followed the 
convention of observational papers, where velocities are de¬ 
fined as being negative for outflow (blue-shifted lines) and pos¬ 
itive for infall (red-shifted lines). 
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4.1. Visual and bolometric phases 

Observed RVs are usually linked to visual phases <py within the 
lightcycle. For the model calculation we know the bolometric 
phases from the luminosity lightcurve (Fig. 01. The question is 
how these are related. For O-rich Miras a small shift is sus¬ 
pected, such that bolometric phases lag behind visual ones. 
Lockwood & Wing (I197H found a typical shift of A(/i«j 0.1 
(if the phases measured by the filter at 1.04jum are taken to 
correspond to bolometric ones). Such a comparison has not 
been published for S Cep or other C-rich Miras yet. In Sect. l4.5l 
a probable shift is discussed from a comparison of line pro¬ 
files and RV-curves with observations. The shift is not ex¬ 
pected to be very large, as the features dominating the V- 
band for C stars (CN, C 2 ) originate in deep photospheric lay¬ 
ers with Tgas =s2500^000K (Fig. 3 of Loidl et al. 1 19991 1 and 
are less temperature sensitive than e.g. TiO in the O-rich case. 
Therefore the V-lightcurve should be closely coupled to the 
pulsation and the bolometric lightcurve. Also, lightcurves of C 
stars are more symmetric and possible phase shifts due to dis¬ 
tortion should be minimal. An additional uncertainty may come 
from the fact that pulsation is simulated by a piston for the dy¬ 
namic model atmospheres and does not come from a pulsation 
model (cf. Hofner et al. I2003t . In the following, bolometric 
phases 0boi will be used to characterise synthetic spectra and 
derived RVs, with numbers written in italic for a clear distinc¬ 
tion from visual phases (py used to denote observations. 

4.2. Probing the pulsational layers 
4.2.1. CO Av=3 lines 

The second overtone CO lines were used most frequently to 
investigate deep photospheric layers in AGB stars and they 
show a very typical behaviour as already described in detail 
in Sect. 5.4.1. of Paper I. Although some observations of these 
lines were reported by Barnbaum & Hinkle J1995II . identifi¬ 
cations are difficult in C star spectra due to contamination by 
other molecules (CN, C 2 ). No time series studies have been 
published yet. The FTS spectra for S Cep of HB96 were all 
of poor quality in the H-band. Therefore, the synthetic line 
profiles in Fig. 4 of Paper I can only be compared to obser¬ 
vations of M- and S-type Miras. The behaviour also is sup¬ 
posed to be similar for C stars. Fig. 1 of Hinkle et al. in 

the following HHR82) shows averaged line profiles for differ¬ 
ent phases for the S-type Mira;^fCyg, which will be used for 
the comparison here. A similar comparison was presented in 
Nowotny et al. J2004t . 

To match profiles of the same shape, a difference in phase 
of «i0.25 between py of the x Cyg observations and 0boi of the 
model is needed. Line strengths are comparable and the syn¬ 
thetic profiles reproduce the typical pattern reasonably well, at 
least for the lower spectral resolution. However, the profiles do 
show some substructures at the spectral resolution of 300 000 
(Paper!). The shifts in wavelength are too small compared with 
observations and the doubling is less pronounced (compare 
phases (f)y-Q.Ql/(pboi-0.75), but still clearly visible. This indi¬ 
cates that the difference in velocity behind and in front of the 



phase 


Fig. 2. Radial velocities for (pboi-0.0-1.0 as derived from syn¬ 
thetic profiles of CO Av=3 lines, shown in Fig. 4 of Paper 1. 
Plotted for spectral resolutions of 300000 (+) and 70000 
(squares) respectively. 


shockfront is too low in the atmospheric model (see Sect. l6.2t . 
A blue-shifted emission feature appearing before the line dou¬ 
bling phases (0.86/0.47) can also be seen in the modelled pro¬ 
files. 

Figure|2l shows RVs derived from these CO lines for var¬ 
ious phases during one pulsational period. The developement 
of a few spectral components can be followed at the higher 
resolution of 300000. A weak component with RVsiOkms ' 
seen at phases of 0.6-0.7 may also (weakly) be recognised in 
the observed Cyg profiles (phase 0.90 in Fig. 1 of HHR82) 
as well as in the composite RV-curve of Lebzelter & Hinkle 
(120021 Fig. 1). It is conspicuous that line splitting is not only 
found around luminosity maximum, but also (but more mod¬ 
erately) during phases ^0.2-0.5, which is not being observed. 
This splitting turns into a continuous transition from blue- to 
red-shifts for a lower resolution of 70000, as the two com¬ 
ponents of the line profile melt into one broad feature (phase 
0.30). Generally, the observed discontinuous RV-curve - com¬ 
pare Fig. 12 of HHR82 - is reproduced by the low-resolution 
profiles. Only the RVs from the lower resolution were thus used 
in the following comparison (Figs. 0and0. The behaviour pat¬ 
tern is highly periodic (Fig.0, as found in observations and 
as also is expected for the model (Paper!). Taking into ac¬ 
count the above-mentioned shift of » 0.25 between (py and (pboh 
the ”S-shape” is reproduced as well as the asymmetry^ w.r.t. 
RV^Okms '. Line doubling appears during phases ^0.7-0.85, 
a similar time interval as observed (A(p^Q.2). Zero-crossing 

^ RV-curves extend to more positive values and appear ’’shifted” to 
RV'=0kms“*. This is seen in all observations of second overtone CO 
lines (e.g. Hinkle et al. 119841 in the following HSH84). As one would 
expect to see equal infall and outflow for radial pulsations, this asym¬ 
metry is interpreted as due to line components that do not come from 
the same depth during the lightcycle. Note that also in the structures 
of the model atmosphere (Fig. 9 in Paper I) infall velocities are larger 
than outflow in the pulsational layers. 
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happens at 0boi=O.i, while it is observed at (p^-QA. Although 
the amplitude^ A/?y»14kms“' is too low (typically 25kms“' 
in Mira observations), the characteristic behaviour could be re¬ 
produced successfully by consistent calculations. 

4.2.2. CN Av=-2 lines 

Figure|3 shows a sequence of FTS spectra of selected phases 
from the most thorough study of S Cep by HB96 with the rel¬ 
atively unblended 1-3 P238.5 line (most other regions of the 
K-band are much more crowded). It reveals that this line shows 
a similar behaviour as CO Av=3 lines. For most of the CN lines 
there, this cannot be recognised clearly due to heavy blending. 
As it was noted by HB96 and was also found for RLeo by 
Hinkle & Barnes M9191 . the IR lines of Ti (e.g. in Fig.|3 at 
4496.625 cm“^) also sample pulsational layers and appear dou¬ 
bled around visual maximum, too. 



Fig. 3. Time series of high-resolution FTS spectra of S Cep 
(from HB96) showing the temporal behaviour of CN Av=-2 
lines of the red system. The rest wavelengths taken from the 
literature are marked. 

In Fig.|3] a comparison between the mentioned observed 
and synthetic spectral lines is presented. Although not ex¬ 
actly the same line - observed CN 1-3 P238.5 vs. synthetic 

^ maximum difference in RP of the two components during phases 
of line doubling 


SCep observed 
R=70.000 



CN A^=-2 synthetic 
R = 70.000 



Fig. 4. Evolution of the observed CN 1-3 P238.5 line (taken 
from the FTS spectra of HB96) compared to synthetic line pro¬ 
files for the CN 1-3 Q24.5 line; both can be found in the K-band 
around 2 yum. 


CN 1-3 Q 24.5 - the behaviour is expected to be very simi¬ 
lar. Only the lower-resolution synthetic spectra from Fig. 5 of 
Paper I are plotted here. The plotted profiles were chosen to 
show maximum congruence, which occurs not necessarily at 
the same phase/phase-shift. As already described in Paper I, the 
change with phase (red/blue-shifts, line doubling) can in princi¬ 
ple be recognised from the synthetic profiles. While the ones at 
higher-resolution show more fine structures than known from 
observations, complexity is reduced by rebinning the spectra 
to lower resolutions and they become more similar to the ob¬ 
served ones. Line strengths are comparable, but again the shifts 
in wavelength are too small. The splitting of the two compo¬ 
nents around maximum phases is much too weak (0.00/0.97). 
They merge and form one broad feature, while observations 
show two distinct lines. This is even more pronounced here 
than for CO Av=3 lines in Sect. l4.2.T1 

Figure|3shows RVs derived from synthetic CN lines during 
the first pulsational period examined. A plot of values measured 
from single profiles at highest spectral resolution looks even 
more complex than for CO lines due to the multi-component 
line profiles. Some components obviously present as a distor¬ 
tion of the line profile do not give a clearly measurable ve¬ 
locity. Opposite to the expectation (motivated by observed and 
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Fig. 6. Synthetic spectrum (resolution R=300 000) dominated 
by CN 0-2 lines of the red system (here for phase 0.0 without 
velocities taken into account), which was used as template for 
deriving radial velocities by cross-correlation. 


Fig. 5. Radial velocities for <p^a\-0.0-l.0 as derived from syn¬ 
thetic prohles of CN Av=-2 lines, shown in Fig. 5 of Paper!. 
Measurements of single lines in spectra with resolutions of 
300000 (+) and 70000 (•) as well as results of the cross¬ 
correlation («) are plotted. 

R=70000 synthetic spectra) of one component (appearing at 
blue-shifts, moving toward red-shifts and disappearing again), 
the development of several components, which appear to be 
independent of each other, can be followed. But profiles re¬ 
binned to a lower resolution of 70000 lead to RV-curves (•) 
that are in notable qualitative agreement with S Cep observa¬ 
tions (Fig.Q. Even the discontinuity in the blue-shifted com¬ 
ponent at 0boi=O.75 is smoothed out. 

Another way to extract velocity curves is not to calcu¬ 
late one single line prohle but rather a whole region occu¬ 
pied by CN lines and derive RVs by cross-correlation. This 
method is not able to resolve very weak components recog¬ 
nisable in single profiles as blended and unblended lines mix in 
the wavelength range used. HB96 used this technique because 
of the crowded CN spectra. We have tried to simulate this ap¬ 
proach with synthetic spectra. Figurel^ shows the spectral re¬ 
gion used, occupied by many CN Av=-2 lines of the red sys¬ 
tem. It was chosen to have not too strong blending, but rather 
several strong individual lines. The positions of CN lines in the 
line list were not corrected for wavelength shifts between theo¬ 
retical and true values; a direct comparison with any observed 
spectra may therefore not be straightforward. Relative veloc¬ 
ities should still be correct. Also for this region spectra with 
resolutions of 300000 were calculated for all phases and the 
line shapes look exactly like the ones in Fig.|4] RVs were then 
measured by cross-correlating these with the template spec¬ 
trum shown in Fig.|3((^boi=0.0, no velocities taken into account 
in the radiative transfer). Results are also plotted in Fig.|3 The 
common characteristic of line doubling can be recognised and 
values similar to the ones from the low-resolution profiles are 


derived for phases (pbo\~0.25-0.8. However, cross-correlation 
is not able to resolve the full amplitude in splitting for phases 
0.8-1.25; the blue-shifted component (RV^-2.5 km s ' instead 
of -7) is especially suppressed. 

As they are more similar to observed ones, only RVs from 
low-resolution single profiles will be considered in the follow¬ 
ing and are plotted in Figs. |Hland0 the latter demonstrating the 
highly periodic behaviour. The observed RV-curve of S Cep in 
Fig. n can be reproduced qualitatively by the synthetic ones. 
The scenario of line doubling due to shockfronts can be fol¬ 
lowed by CN line profiles, too. No global shift (as for CO 
Av=3 lines) between observed and (^boi of the models is 
needed to match the RV-curve of S Cep and the model. Line 
doubling appears almost for the same interval (0boi~ 0.85-1.1). 
Zero-crossing occurs at the same phase (0.5). The synthetic RV- 
curve is symmetric w.r.t. RV^Okms''. The velocity amplitude 
of ARyssl3.5kms ' is too low compared with the observed 
one of 22 km s ' for S Cep. Similar large amplitudes from CN 
lines in the optical have also been reported for a few other 
Miras by Barnbaum J1992t . 

4.3. Probing the dust-forming region - CO Av=2: 

Figure0 shows a comparison of first overtone CO lines of 
S Cep observations (HB96) with synthetic profiles, already pre¬ 
sented and described in detail in Paper I. The weaker high- 
excitation lines (e.g. R82) exhibit the same behaviour as CO 
Av=3 lines in the synthetic spectra, which can be difficult to 
see in observed C star spectra due to blending with several 
other lines. The very strong low-excitation lines (R18,19) can 
be identihed clearly though and are studied in more detail here. 
While the prohles differ somewhat for different lines (compare 
e.g. R19 and R18) in spectra observed at a given time (due 
to contamination and observational uncertainties), they are ex¬ 
actly reproduced for synthetic ones. Line depths are compa¬ 
rable for synthetic and observed spectra as well as the general 
appearance (see also Lig. 3 in HB96). These lines have complex 
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Fig. 7. Upper panel: Observed CO Av=2 lines from the FTS 
spectra of S Cep (from HB96) for selected phases. 

Lower panel: Synthetic line prohles (f?=300000) for compar¬ 
ison. The low-excitation lines (2-0 R19 and R18) consist of 
several components and appear always blue-shifted, while the 
high-excitation line (2-0 R82) follows the periodic behaviour 
of CO Av=3 lines. 


shapes with various components present, although no clear line 
splitting is seen most of the time. Line prohles change with 
time, but these changes are not coupled to the lightcycle. This 
complex behaviour has also been found in obervations of R Leo 
by Hinkle lfT978t and of Cyg by HHR82 (their Fig. 8). Prohle 
variations in the CO hrst overtone low-excitation lines that are 
not related to the light-cycle were also found earlier, e.g. by 
Winters et al. J2000at . who used a very similar approach. 

Radial velocities were derived from the deepest points of 
the synthetic line prohles (at a resolution of 300 000) and are 
shown in Figs.and(open circles). Weak line doubling was 
found only in some cases, disappearing in rebinning to resolu¬ 
tions of 70000 (Fig. 7 in Paper!). Sampling the region of the 
onset of the wind, the main component of low-excitation CO 
Av=2 lines appears slightly blue-shifted by »5kms“* most of 
the time. In observed spectra, RVs for these lines were found 
to be either roughly at the CMRV ix Cyg) or show some vari¬ 
ation (RLeo, S Cep or Lebzelter et al. ll999t . but with smaller 
amplitudes than lines sampling pulsational layers (CO Av=3, 
CN, CO Av=2 high-excitation). These three behaviours can be 


understood from the velocity structures of Fig. 9 in Paper!, if 
the lines probe layers of 3, 2 and 1.5 R* respectively. For a!! 
stars studied, RVs from CO Av=2 lines were not coupled to 
the lightcycle. This may rehect the irregular motions within the 
dust-forming region (Sect. 4 in Paper I), where CO Av=2 lines 
originate. 

4 . 4 . Probing the outflow - CO Av=1: 

As described in Paper!, lines of fundamental CO bands are us¬ 
able to study the outflow of mass-losing Miras. Unfortunately, 
all of the few existing LM-band spectra for S Cep are of poor 
quality. A comparison can therefore only be done with other 
results. HHR82 presented in their Fig. 11a spectrum of x Cyg 
showing strong fundamental CO lines at 4.6pm. Being affected 
by several photosperic lines as well as telluric CO lines, they 
still show a constant blue-shift over time. The same lines were 
observed in spectra of IRCH-10216 by Keady et al. and 

later compared to models by Winters et al. OOOOat . Prohles 
there compare well to the ones shown in Fig. 8 of Paper I, 
except for the saturated absorption. They show PCygni-type 
shapes, which is typical^ for lines sampling the outermost re¬ 
gions and indicating stellar winds. RVs were calculated from 
the minima of the blue-shifted absorption in prohles with high¬ 
est resolution (300000); the results are shown in Fig.|3 The 
sometimes very broad or asymmetric absorption features can 
complicate the RV measurements. Apart from some small long¬ 
term variation in the terminal velocity, CO Av=l line velocities 
show a steady outhow over time. 

4 . 5 . The overall picture 

Figurel^shows a compilation of all RVs derived from synthetic 
line prohles presented above. Calculations were done for three 
separated periods (23 instances of time each). Figure|3 shows 
a plot for direct comparison with the observations of S Cep in 
Fig-ID For CO Av=3 and CN Av=-2 lines RVs from all pe¬ 
riods were combined into one composite light-cycle and then 
plotted repeatedly for better illustration (RV-curves from differ¬ 
ent periods are almost identical each other, because of the very 
regular movement of the inner regions, see Fig. 2 in Paper I). 
For CO Av=2 and Av= 1 lines only measurements from one pe¬ 
riod (0boi=7-S and 14-15 respectively) are plotted for a clearer 
picture. It can be recognised that (with restrictions) our model 
calculations can reproduce the velocity pattern found for S Cep 
(and other Miras). 

The hrst two types of lines sample deep pulsating layers 
and show typical discontinuous RV-curves, rehecting shock- 
fronts running through the line-forming regions. Synthetic 
line prohles could generally reproduce this characteristic be¬ 
haviour, at least for the lower spectral resolution of 70000 
(only these results are adopted here). Some aspects (S-shape, 

* This is however not necessarily the case for all lines being formed 
in the outflow. CO Av=2 lines for example, although sampling the re¬ 
gion where the wind starts, are not expected to show P Cygni-emission 
as the de-excitation v=2 —»0 has a lower transition probability than 
cascading via the v=l-state. 
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Fig. 8. Upper panels: Lightcurve of the dynamic model atmosphere for three different, separated periods. The phases <^boi for 
which synthetic spectra were calculated and the reference phase 0.00 are marked. 

Lower panels: Compilation of the corresponding radial velocities derived from shifted line prohles for different types of molecular 
lines. Note the velocity convention: Negative velocities (blue-shift) represent outflowing matter, positive values (red-shift) denote 
infall. 


asymmetry w.r.t. RV—0, line doubling interval) are realistically 
replicated with synthetic CO Av=3 lines, but the apparant phase 
shift is conspicuous. While the mean phase of the line doubling 
interval is 0v~l.O in most observations (HSH84), it is clearly 
shifted to (p\,oi~0.775 in Fig.|3 suggesting a shift of A0ssO.225 
between visual and bolometric phases for our models, such that 
0boi lags behind 0v- Such a phase shift is not found for the RV- 
curve of CN lines however, which resembles the observations 
of S Cep in Fig.^ rather well. Line doubling appears for the 
same interval and at the same phases, as well as zero-crossing. 
The synthetic RV-curve is symmetric around RV-Q (in contrast 
to CO lines), but the observed one even extends to more neg¬ 
ative values. CN lines appear to originate slightly further out^ 


^ Combining the difference of the median phases of line doubling 
for both type of lines (A0boi~O.975-O.775) with an estimation for the 
shock front propagation velocity from Fig. 2 of Paperl (dR/d0boi~O.88 


than the CO second overtone lines, which can be inferred from 
plots of the radial gradient of t (as Fig. 3 in Paper I) or from 
line doubling at later phases (the shockfront propagating out¬ 
wards reaches outer layers later, also demonstrated in Fig. 4 of 
Alvarez et al. I2000t . Such a comparison would be interesting 
if done observationally in the future. Dehnitely, the amplitudes 
ARV of both types of lines are too small compared to observa¬ 
tions (see Sect. l6.2> . 

Although synthetic prohles of CO Av=2 low-excitation 
lines are comparable to observed ones, absolute values of RVs 
are qualitatively different. They show a small blue-shift at all 
times and therefore sample the region where the stellar wind 
starts («i3 R*). On the other hand, observed RVs are either con¬ 
stantly at siCMRV or vary little around the CMRV. Comparing 

^ Mfront~7.2kms“') this would correspond to a difference AR of 
0.176R./86.8Ro/0.4AU. 
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Fig. 9. Radial velocities from Fig.|8l projected onto one light- 
cycle. While CN and CO Av=3 lines probe deep pulsating lay¬ 
ers, CO Av=2 lines come from the dust-forming region where 
the stellar wind is triggered and the outflow is sampled by 
CO Av=l lines. To be compared with observational results for 
S Cep in Fig.^ 

Fig. 9 of Paper I, this could be interpreted as either the lines 
should probe deeper layers (a;2 or 1.5 R*) or pulsation should 
influence layers at larger radii. 

Sampling the steady outflow, CO Av=l lines appear clearly 
blue-shifted, more or less constantly over time. Compared to 
observed RVs from KI lines in S Cep spectra (Fig.[0, the syn¬ 
thetic ones are ss5 km s ' too low (as already expected from the 
velocity structures of the model, see Fig. 9 in Paper I). 

The velocity stratigraphy sketched for x Cyg around max¬ 
imum light in Wallerstein (fT^ can directly be recognised 
in the movement of mass shells of the model atmosphere for 
<Pho\~0-0. 

4.6. Gas velocities vs. measured RVs 

To infer actual gas velocities from RVs of shifted spectral fea¬ 
tures, a correction factor p - with Ug^^-p-RV - should be con¬ 
sidered, caused by geometric aspects and center-to-limb ef¬ 
fects. A first theoretical estimation of 24/17=sl.41 for any ra¬ 
dial pulsating star was given by Getting (|I935}. Parsons 
later found smaller values from synthetic lines calculated with 
Cepheid model atmospheres and a dependence on chosen pa¬ 
rameters (his Fig. 2). His approximate value of 1.31 was then 
used by Hinkle to analyse the observations of RLeo. 

As the very extended Mira atmospheres have more compli¬ 


cated velocity structures (e.g. Fig. HU) and lines from different 
wavelengths sample different depths, this correction becomes 
more difficult. As stated in Scholz & Wood | |2000> a factor of 
ssl.4 was often used to interpret the observations. Comparing 
RVs derived from line shifts with gas velocities in front of and 
behind the shockfronts of their models, they derive factors of 
1.25/1.4/1.6 for different lines (CO, OH) from model calcula¬ 
tions. Using the same method, a correction factor of p » 1.36 
from synthetic CO Av=3 lines at phases 0.75/0.85 (Sect. l4.2Tl 
could be estimated here. 

5. Quasi-static, warm moiecuiar enveiopes and 
dynamic modei atmospheres 

Different velocity components have been found by studies of 
line profiles in observed high-resolution spectra of pulsating 
AGB stars in the past (HHR82). Among these, a layer with al¬ 
most zero velocity relative to the stellar systemic velocity and 
an excitation temperature of about 800- 1000 K could be iden¬ 
tified from low-excitation lines of CO (Av=2) and OH. Being 
present in a number of Miras [e.g. in RLeo (Hinkle fl978t or 
xCyg (HHR82)], this layer is variable in intensity. The vari¬ 
ations do not follow the lightcycle and sometimes the static 
layer vanishes completely. HHR82 interpreted this layer as a 
reservoir for mass loss driven by radiation pressure on the dust, 
which can condense at these low temperatures. 

Tsuji ( I1988> found anomalies in low-excitation CO Av=2 
lines in high-resolution spectra of various M giants (with 
slightly higher temperatures and smaller light amplitudes). 
These could not be explained by photospheric absorption alone 
but seemed to be composed of two components. Subtracting 
synthetic spectra (based on hydrostatic model atmospheres) 
from observed ones revealed excess absorption with small 
blue- or red-shifts. Tsuji also interpreted this as being caused 
by a ’’quasi-static layer” in the outer atmosphere and derived 
excitation temperatures of 1000-2000 K as well as mass es¬ 
timates for it. Such a molecular formation zone could be well 
separated from the warm photosphere but not yet in the expand¬ 
ing regions of the cool stellar wind. Applying the same method 
to ISO SWS spectra, Tsuji et al. J1997> found excess absorp¬ 
tion/emission for other spectral features (molecular bands of 
H 2 O, CO, SiO, CO 2 ), too. They claimed that such a ’’warm 
molecular layer” could be located at about 2 R, and could be 
common not only in Miras, but in evolved red giants in gen¬ 
eral. Still it was not clear how such an extra component within 
a red giant atmosphere could develop and persist. 

The atmospheric model used here shows a similar be¬ 
haviour to that described above (compare Figs. 2 and 9 in 
Paper I). Around 2 R, we find a region where gas velocities are 
rather small at all times. In this transition zone between the pul¬ 
sating photosphere and the cool outflow, dust is being formed 
(cf. Fig. Ic in Paper!) and the stellar wind is triggered. Low- 
excitation CO Av=2 lines, which sample these layers, show 
only small blue-shifts and RVs would be observable close to 
the stellar systemic velocity. This means that dynamic model 
atmospheres of Hofner et al. J2003t can produce a phenomenon 
resembling a quasi-static layer in a consistent way, and it is 
not necessary to introduce it artificially. This is at least true for 






























10 


Nowotny et al.: Atmospheric Dynamics in Carbon-rich Miras II. 


some combinations of stellar parameters. As it will be shown 
in Sect. l6.ll a relatively small change in the chosen parameters 
(e.g. 1^1.5 Mq) can strongly influence velocities in the men¬ 
tioned layers, due to the non-linear response of dust formation 
to the thermodynamic conditions in the upper atmosphere. 

It will be interesting to perform a parameter study of var¬ 
ious models with regard to this phenomenon. On the other 
hand observational evidence of a static layer for a statistically 
meaningful number of stars should be possible with the up¬ 
coming IR spectrographs. There is already some effort being 
spent to investigate existing high-resolution spectroscopic data 
(Lebzelter et al. In addition, interferometric measure¬ 

ments of Miras will help to tackle this question (e.g. Perrin et 
al.BICT. 

The possible formation mechanism of ’’warm molecu¬ 
lar layers” was investigated by Woitke et al. lfT9^ for O- 
rich Miras and by Helling & Winters ll2()()lal 1200 lb> for 
C-rich ones. They found enhanced molecular column den¬ 
sities in the layered structures of dynamic models of AGB 
outer atmospheres (compared with predictions of hydrostatic 
models), caused by stellar pulsation and propagating shock 
waves. Similar results regarding the partial pressures of se¬ 
lected molecules were described by Hron et al. dT^ and 
Loidl et al. ([T 999 J, also discussing the effects on the resulting 
synthetic spectra of C-rich models. 

6. Future intentions 

6.1. Fitting modeis to observations 

As pointed out in Sect. l4.5l our model can qualitatively repro¬ 
duce the velocities found for S Cep, but some differences re¬ 
main concerning the absolute values. Therefore we tried to vary 
the input parameters of the model (Table 1 of Paper I) to get 
even more realistic atmospheric structures. 

For example a model where only stellar mass is changed 
from 1 to 1.5 Mq was investigated. This seems comparably 
plausible, as we only have rough mass estimates for S Cep (e.g. 
2.5^Mo from Barnbaum et al. I1991> . Differences are imme¬ 
diately visible from the structures plotted in Fig. cni (compare 
Figs. 1 and 9 in Paper I). Dust shells and density enhancements 
are more strongly pronounced. The global velocity structure 
looks markedly different. The very regular motions due to pul¬ 
sation only reach out to ss 1.3 R*. From there on the model struc¬ 
tures are not necessarily similar for the same phases of different 
lightcycles any longer, because the dust formation cycle spans 
two pulsational periods. A plot with movements of mass shells 
(as Fig. 2 in Paper I) looks similar to Fig. 2b in Hofner et al. 
( I2003> . Throughout the dust-forming region strong variations 
in velocities can be seen; a steady outflow and a smooth distri¬ 
bution of the degree of dust condensation is not reached until 
^lOR.. 

These differences are also reflected in line profiles and RVs. 
While the RV-curve for CO Av=3 lines is almost the same 
as for the original model, line strengths change considerably 
from one period to the next (duplicating only every second pe¬ 
riod). CN Av=-2 lines are more affected by the double peri¬ 
odicity. Not only the line shapes (variable line strengths; line 


doubling much clearer), but also RV-curves (that are irregu¬ 
lar with a larger ARV of sil6kms“' for cross-correlation) evi¬ 
dently deviate. However, synthetic velocity amplitudes for lines 
from the inner regions could not be increased enough with this 
changed model to reproduce observed ones. While the faster 
outflow leads to more realistic RVs measured from CO Av=l 
lines («;18kms“'), the largest differences by far compared to 
the original model are found for CO Av=2 low-excitation lines. 
As suspected from the large variations of gas velocities within 
the dust-forming region at *2-3 R* in Fig. cni these lines show 
various components, in part with even more pronounced vari¬ 
ability than second overtone lines (which has never been ob¬ 
served so far). 

From a few other studied models with slightly changed 
characteristics, it appears that the fine-tuning to achieve a ded¬ 
icated fit for the velocities of one particular star (like S Cep) 
is not a straightforward task (compare also the discussion in 
Sect.|3- Not only are the stellar parameters of S Cep poorly 
known (L, M, etc.), but a multidimensional grid of dynamic 
models would also be needed. Too many parameters interact 
and influence the velocity structures in subtle ways. Especially 
the transition region between regular pulsation and steady out¬ 
flow (and thus the behaviour of first overtone low-excitation 
CO lines) seems to depend sensitively on the chosen set of 
stellar parameters. This dependency was discussed also by e.g. 
Hofner & Dorfi ^\WT\ or Winters et al. J2000bt . 


6.2. Larger veiocity ampiitudes 

The initial aim of future calculations will be a dynamic model 
atmosphere with a more realistic velocity amplitude within the 
pulsational layers. As stated in Sect. 2.2 of Paper I, the RV- 
curves of CO Av=3 lines show a rather uniform picture for all 
Miras studied so far, independent of their properties. Figure 1 
of Lebzelter & Hinkle J2002t and the results in HSH84 demon¬ 
strate that amplitudes of ARyai20-30kms“* can be observed. 
This common feature seems to be a fundamental characteris¬ 
tic of Miras and should be reproduced by a realistic model. 
The ARV values given in Scholz & Wood (I2000t could serve 
as a guideline for subsequent modelling. They state that a typ¬ 
ical Mira appears to have a true (observed RVs corrected by 
their factor p to get real gas velocities, Sect. l4.6> post-shock 
outflow velocity of «;-14kms ' and a pre-shock infall ve¬ 
locity of a!20kms ', resulting in a full velocity amplitude of 
a:34kms '. Comparing this with gas velocities in the inner re¬ 
gions of the model used here (Fig. 9 in Paper!), it can be de¬ 
duced that the dynamic model should have an amplitude larger 
by a factor of ^2. With these more extreme velocity gradi¬ 
ents across the shockfront, line doubling should be more pro¬ 
nounced for CO Av=3 and CN lines and easier to recognise 
even in low-resolution spectra (Sect. O- From our studies so 
far we know that a larger shock amplitude cannot be achieved 
simply by applying a piston with a higher velocity amplitude 
Aup. This is due to the self-regulating interdependence of the 
density structure (levitation) and the maximum shock strength. 
More detailed work is needed to address this problem. 
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Fig. 10. Radial structures of a dynamic model atmosphere with 
parameters as in Table 1 of Paper I, except for the stellar mass, 
which is changed to 1.5 Mq. To be compared with structures of 
the original model in Fig. 1 of Paper I. 


6.3. 0-rich model atmospheres 

Calculations presented here and in Paper I focus on carbon-rich 
model atmospheres, as the formation and evolution of C-rich 
dust is much better understood and implemented in these mod¬ 
els. More data have been obtained observationally for M-type 
stars. Thus further investigations of synthetic line profiles (for 


CO lines but also for others found in spectra of O-rich stars, as 
OH, SiO, etc.) based on oxygen-rich models may be revealing. 

7. Conclusions 

From high-resolution spectroscopic studies it is known that dif¬ 
ferent molecular absorption lines in the NIR sample different 
depths of the dynamic atmospheres of pulsating AGB stars. 
Gas velocities there can heavily influence line profiles (Doppler 
shift, broadening, line doubling). Reproducing this scenario 
by consistent calculations is an indicator of the quality of dy¬ 
namic model atmospheres. The advanced models of Hofner et 
al. J2003t allow us to simulate the dynamics in C-rich Mira 
stars from the innermost layers dominated by pulsation out to 
the regions of the stellar wind. Synthetic line profiles (pre¬ 
sented in Paper I) were compared in detail with observations 
of S Cep and other Miras. In addition, radial velocities were 
derived from their shifts in wavelength. 

It was found that the characteristic behaviour of spectral 
features sampling different regions (pulsational layers, dust¬ 
forming region, outflowing outer atmosphere) could qualita¬ 
tively be reproduced by consistent calculations. Line profiles 
compare reasonably well to observed ones. Photospheric lines 
of CO (Av=3) and CN (Av=-2) appear more complex at high¬ 
est spectral resolution. Rebinning to lower resolutions (com¬ 
parable to observed FTS spectra) leads to the expected, very 
typical S-shaped RV-curve. Together with RVs of low excita¬ 
tion CO Av=2 lines (sampling layers where the stellar wind 
is triggered) and of CO Av= 1 lines (probing the outermost re¬ 
gions with strong outflow) we made a direct comparison with 
the equivalent data for S Cep. This comparison reveals that the 
dynamic model atmosphere used can reproduce - at least qual¬ 
itatively - the global velocity structure throughout the atmo¬ 
sphere of a C-rich Mira. 

Nevertheless, there are still quantitative differences: for ex¬ 
ample both the RV amplitudes in the inner regions and the out¬ 
flow velocities are too small. Generating a model with quan¬ 
titatively fitting radial velocities for a particular observed star 
will be challenging (Sect. l6.ll . A realistic model should have a 
larger velocity amplitude in the pulsating photosphere, which 
appears to be a common feature of Mira stars (Sect. l6.2l . 

The interplay between pulsation and the onset of a dust- 
driven wind can result in a transition zone with rather low gas 
velocities, at least for models with certain parameters. Thus, 
the dynamic model atmospheres used can provide a natural 
and phyically consistent explanation for a supposed quasi-static 
layer within the atmosphere, which was proposed for some red 
giants from spectroscopic observations (Sect.|3. 

In addition to high-resolution spectroscopy, interferometric 
(narrow band) observations with high spatial resolution (e.g. 
Hron et al. 120031 will be a complementary major tool for stud¬ 
ies of Mira atmospheres in the future. Time series of syn¬ 
chronous observations with both methods could provide un¬ 
precedented insights and strong constraints for atmospheric 
modelling. 
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